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ABSTRACT 

New single dish and millimeter wave interferometer observations of the molecular ab- 
sorption line system in the gravitational lens to PKS1830-211 at z=0. 88582 is presented. 
Self-calibrated interferometer data shows unequivocally that the previously detected ab- 
sorption component is associated with the gravitationally lensed south-west image of the 
background source. A second absorption line of HC0"*'(2^1) at z=0. 88582 is detected. 
This component is shifted in velocity by -147kms~^ relative to the main absorption 
line, and is shown to be associated with the north-east image. A few additional ab- 
sorption lines are presented. Upper limits to absorption and emission lines from the 
possible absorption system at z=0.1927, seen in 21cm HI by Lovell et al, are reported. 
Alternatively, there could be only one absorbing system at z=0. 88582, with the HI line 
attributed instead to the PS molecule. 

Subject headings: ISM : molecules — BL Lacertae objects: individual (PKS1830-211) 
— galaxies : abundances — galaxies : ISM — quasars : absorption lines — gravitational 
lensing 
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1. Introduction 

The molecular interstellar medium (ISM) in galax- 
ies at rcdshifts up to z»0.9 has been studied through 
the use of absorption of molecular rotational tran- 
sitions (cf. Wiklind & Combes 1996a, b; Combes & 
Wiklind 1995). Detection of several different lines 
of molecules such as CO, HCO+, HCN, HNC, CS, 
N2H+, H2CO, etc. allows a derivation of accu- 
rate excitation temperatures, column densities and 
abundance ratios (Wiklind & Combes 1994, 1995, 
1996a,b). Four molecular absorption line systems are 
known to date; in two of them the absorption occurs 
in the host galaxy to the AGN and two represent truly 
intervening galaxies. Since the distribution of molec- 
ular gas in galaxies is usually centrally concentrated, 
the likelihood to encounter a molecular cloud is higher 
if the line of sight passes close to the center of an in- 
tervening galaxy. This means that also the likelihood 
for gravitational lensing is high and, indeed, in each of 
the two systciuis where absorption occurs in an inter- 
vening galaxy, that galaxy acts as a gravitational lens 
to the background source (Wiklind & Combes 1995, 
1996a) 

The use of molecular absorption lines is manifold. 
Apart from studying the physical and chemical prop- 
erties of molecular gas in distant galaxies, they can 
be used as probes of very small scale structures in 
the molecular gas. They can also be used as cosmo- 
graphic probes, to set upper limits to the temperature 
of the cosmic microwave backgroimd radiation and to 
constrain the geometry in gravitationally lensed sys- 
tems. The molecular ISM is probed on a scale of a 
few tenths of a parsec, limited only by the angular ex- 
tent of the background source. The main uncertainty 
in the molecular gas properties derived from the ab- 
sorption lines is due to our ignorance of possible small 
scale structures in the ISM on scales smaller than the 
angular extent of the background source (e.g. Wik- 
lind & Combes 1997). 

In this article, we concentrate on the absorption 
line system originating in a galaxy which plays the 
role of a gravitational lens at z=0.886 (Wiklind & 
Combes 1996a). The galaxy absorbs the continuum of 
the strong flat-spectrum radio source PKS1830-211. 
This radio source has long been considered a gravita- 
tional lens candidate on basis of the radio structure 
(Rao & Subrahmanyan 1988; Jauncey et al. 1991). 
Since optical and near infrared searches for an optical 
counterpart to the radio source as well to the lensing 



galaxy have remained unsuccessful (cf. Djorgovski et 
al. 1992), the modeling of the lensing system was lack- 
ing redshift information (Kochanek & Narayan 1992; 
Nair, Narasimha & Rao 1993). The galaxy at z=0.886 
was discovered by searching for molecular absorption 
lines in the 3- and 2-mm band (Wiklind & Combes 
1996a). A second intervening system has been re- 
ported at z=0.193, seen through 21 cm HI absorption 
by Lovell et al. (1996). 

We present new molecular absorption data on the 
galaxy at z=0.886 as well as upper limits to molecular 
absorption and emission on the z=0.193 system. We 
also discuss what can be inferred regarding the mass 
of the lensing galaxy and the redshift of the back- 
ground source. 

2. Observations and data reduction 

Our data were obtained using the IRAM millime- 
ter wave interferometer situated on Plateau de Bure 
in Prance, the single dish IRAM 30-m telescope on 
Pico Veleta in Spain and the 15 m SEST (Swedish 
ESO Submillimeter Telescope) on La Silla in Chile. 
The interferometer observations were done on several 
occasions during 1995 and 1996. The single dish ob- 
servations were likewise obtained during several ob- 
serving runs in 1995 and 1996. 

At the IRAM 30-m telescope we used 3-, 2-, and 
1.3-mm SIS receivers tuned to the redshifted fre- 
quency of the observed transition. The observations 
were made using a nutating subreflector with a fre- 
quency of 0.5 Hz and with a throw of ±90" in az- 
imuth. Typical system temperatures in the scale 
were 400 600 K. The weather conditions were good 
on most observing sessions. Pointing, focusing and 
calibration were done before each observation. As 
backends we used narrow band autocorrelators, giv- 
ing a velocity resolution of approximately 0.7kms~^ 
and a bandwidth of 400kms~^ At the SEST we used 
3- , and 2 mm SIS receivers with typical system tem- 
peratures 200-300 K. Here we used an AOS with a 
bandwidth of 1200 km s~^ and a velocity resolution 
of 0.7 km s""'^. These observations were done in a dual 
beamswitch mode, where the beam is displaced 12' 
in azimuth at a frequency of 6 Hz. The SEST data 
was obtained under good weather conditions. Point- 
ing and calibration where done as at IRAM. 

The interferometer data were made with the stan- 
dard BC configuration (see Guilloteau et al. 1992). 
The array comprised 4 15-m telescopes. The receivers 
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were 3-mni SIS, giving a typical system temperature 
of 400-500K. PKS1830-211 itself was used as a phase 
reference, while bandpass and amplitude calibrations 
were done using also NRAO 530, 3C454.3 and W30H. 
The data reduction was made with the standard CLIC 
software. Due to the low elevation of the PKS1830- 
211, we obtained a beam of 6 x 2.4 arcsec. This ap- 
pears insufficient to resolve two images separated by 
only one arcsec; but due to the high signal-to-noise 
ratio, we were able to fit the integrated line profile 
and the phase center, as a function of frequency, di- 
rectly from the visibilities. We also made a two point 
sources fit, which gave the relative positions of the 
two images and their relative fluxes. 

The accuracy of the measurement depends on the 
inherent uncertainty in the chopper-wheel calibration 
method as well as the pointing of the telescope. The 
overall accuracy for a single observation on the inten- 
sities are around 5-10% for the interferometer data 
and 10-20% for the single dish data. 

3. Results 

3.1. The location of the z=0.886 absorption 

In our first study of the molecular absorption 
line system towards PKS 1830-2 11 we found that the 
HC0+(2^1), HCN(2^1) and, possibly, the HNC(2^1) 
lines were heavily saturated (Wiklind & Combes 1996a). 
Yet, the depths of these absorption lines were only 
36% of the total continuum. Since the magnifica- 
tion ratio derived from long wavelength interferom- 
eter data is 1.8 (Nair et al. 1993), corresponding to a 
flux fraction of ~36% for the SW and 64% for the NE 
corej^, respectively, we concluded that the molecular 
gas only covers the SW core, with a high covering 
factor (close to 100%). The present interferometric 
results confirm that only the SW component is ab- 
sorbed around V = 0. 

In Fig. 1 we show the HC0+(2^1) spectrum ob- 
tained with the Plateau de Bure interferometer. The 
low declination of PKS1830-211 makes imaging diffi- 
cult, partly due to a north-south elongated beam, and 
to the large air-mass. The continuum is, however, 
strong enough to allow a self-calibration of the data. 
We made a gaussian fit of the source directly from the 
visibility data, and traced the location of the phase- 
center as a function of frequency. The accuracy is 

Henceforth we refer the the north-east component as NE and 
the south-west component as SW. 



very good in RA (^0.05") and slightly worse in DEC 
(~0.1"). The two top panels in Fig. 1 show the rela- 
tive phase in RA and DEC as a function of velocity 
across the absorption profile. From the data obtained 
in September and October 1995, the maximum shifts 
are Aa = 0.22" and AS = 0.28". From the April 
1996 data we get Aa = 0.24" and AS = 0.26" (see 
also Fig. 3). 

If the flux fraction from the NE and SW cores are 
/i and /2, respectively (/i + /2 = 1), the shift of the 
phase-center along the line joining the cores is 

Aa = (3{v)f2r sintp , AS = (3{v)f2r costfi , 

where f3{v) is the fraction of continuum flux from the 
SW core that is covered by optically thick molec- 
ular gas of velocity v, r is the separation of the 
two components and ■0 is the position angle of the 
line joining the two cores, measured from north to 
east. From the 15 GHz VLA continuum map pre- 
sented by Subrahmanyan et al. (1990), we derive 
r = 0.98" and t/j = 44°. Expressing the flux frac- 
tions in terms of the magnification ratio R — /1//2, 
we get f3 = {1 + R)Aa/{rsmtp). Using the mag- 
nification ratio R derived from long wavelength in- 
terferometry of 1.8 (Nair et al. 1993) and assuming 
that the covering factor /3 is 100% for the SW compo- 
nent, the phase-center should shift Aa « 0.24" and 
AS « 0.25", which is consistent with the shift that 
we see in Fig. 1. 

Frye, Welch & Broadhurst (1997) observed PKS1830- 
211 with the BIMA interferometer and derived a mag- 
nification ratio of 1.14 ± 0.05 by resolving the two 
continuum components. Since the HC0+(2<— 1) line 
is manifestly optically thick, as seen through the de- 
tection of the H"C0+(2^1) line (Wiklind & Combes 
1996a), such a low magnification ratio indicates that 
the covering factor f3 of the SW core is only 80- 
84%. However, they measure a total flux of 2.25 Jy in 
March 96 (1 March and 20 March), much larger than 
ours of 1.55±0.06 Jy on 28 February and 1.45±0.08 Jy 
on 14 April 96). Measurements made with the Owens 
Valley interferometer on 1 March gave 1.6 ± 0.1 Jy 
(S. Aalto priv. comm.). It thus appears that the con- 
tinuum is overestimated by a factor 1.4-1.5 in Frye 
et al. (1997). We stress that calibration is quite deli- 
cate since there exist significant rapid variations in the 
continuum levels of the two images. Given the high 
signal-to-noise of the Plateau de Bure observations, 
we have tried to retrieve the magnification ratio from 
a two point sources fit of the visibility data. The fits 
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were very robust, completely independent of the ini- 
tial parameters, and revealed an RA and DEC sepa- 
ration between the two images exactly consistent with 
r ~ 0.98" and tp = 44° in all the 136 channels with- 
out absorption. In Table 1 we show the magnification 
ratios obtained at the various epochs. Although the 
absolute calibration is not so precise, we believe that 
the relative ratios between the two components is ac- 
curate. Although the magnification ratio has changed 
from 1.8 ±0.1 to 1.2 ±0.1 between October 1995 and 
April 1996, in concordance with the value derived by 
Prye et al, the lower continuum level measured by us 
is in agreement with a covering factor ^1 for the SW 
component (see Sect. 4.2). 

3.2. A second absorbing cloud at z=0. 88489 

We have detected a second absorption component 
situated at a velocity — 147kms~^ relative to the 
main HC0+(2<-l) absorption at z=0.88582. The sec- 
ond component is hereafter refered to as the 'satel- 
lite'. The HC0+(2<— 1) line has been observed ex- 
tensively during a monitoring campaign at both the 
IRAM 30 m and SEST 15-m telescope. Averages of 
data obtained during 1996 are shown in Fig. 2, where 
the two dashed vertical lines mark zero velocity and 
— 147 km s"'^ in a heliocentric velocity scale in the 
restframe of the z=0. 88582 absorber. The satellite 
line is not seen in the HC0+(3<— 2) spectrum. Prye 
et al (1997) could have detected this component in 
March 96 (since the line was present at this epoch, 
see Fig. 3), but their lower signal- to noise ratio pre- 
vented the detection. 

Could the second absorption line be associated 
with the z=0.193 HI absorber rather than the one 
at z=0.886? The only line at z=0.193 which could 
be in the 'vicinity' of the observed frequency of the 
satellite line at 94.635 GHz is the redshifted C0(1- 
0) line at 96.650 GHz. The difference is in excess 
of 2 GHz, corresponding to a velocity difference of 
^6200 km s~^ in the restframe of the z=0.193 ab- 
sorber. It is therefore highly unlikely that the satellite 
absorption is associated with this system. We con- 
clude that the satellite line is a second HC0+(2<— 1) 
line at a blucshiftcd velocity of — 147kms~^ relative 
to the redshift z=0. 88582 (corresponding to a redshift 
of z= 0.884896 ±0.00001). 

Whereas the main absorption line at y = has a 
depth of 38-40% of the total continuum, the satellite 
line has a depth of only 7%. Measured in antenna 
temperature, this corresponds to 1.5 mK and 6mK in 



the SEST and IRAM data, as of April 96. The actual 
opacity, however, depends on the surface covering fac- 
tor of the absorption. The separation of the NE and 
SW contimmm components is 0.98", corresponding 
to a distance of 5.6 kpc, measured in the plane of the 
sky (Ho = 75kms~^Mpc~\ qo=0.5). The satellite 
absorption can therefore be associated with only one 
of the continuum sources. Since the fiux contribu- 
tion from the NE and SW components are ~60% and 
^40%, respectively, the satellite absorption line is ei- 
ther optically thin or has a very small covering factor. 

The satellite absorption line is clearly seen in the 
interferometer data (Fig. 3). The relative location of 
the phase-center for Aa is shown in the top panel. 
In addition to the positive shift of the phase-center 
already established for the main absorption around 
V = 0kms~^, a shift is also seen at the velocity of 
the satellite line, but this time towards negative val- 
ues. This shows that the phase-center moves partly 
towards the SW continuum component at the veloc- 
ity of the satellite absorption. The second absorption 
component must therefore originate in front of the NE 
continuum core. We are thus seeing molecular gas 
from two different parts of the intervening galaxy. 

Two absorption components, separated by 330 km s" 

are also seen at z=0.673 towards the radio source 
B3 1504±377 (Wikhnd & Combes 1996b). In this case 
the absorption occurs in the galaxy hosting the radio 
source. The galaxy is close to edge-on and the sepa- 
ration of the two absorption lines are caused by non- 
circular velocities along a single line of sight within 
the disk of the host galaxy, most likely caused by the 
presence of a barred potential. The line of sight to 
the two images of the continuum core of PKS1830- 
211 passes through the intervening galaxy at ^1.8 and 
~3.8kpc (SW and NE components) from the lens cen- 
ter (Nair et al. 1993). A normal rotation curve (and 
circular orbits) are sufficient here to account for the 
observations. 

We derived column densities using the formalism 
presented in Wikhnd & Combes (1995, 1997). The 
HCO+ column density of the satellite absorber, de- 
rived using 60% of the total flux, is only 7 x 10^^ cm~^ 
(however we can not cxchidc the case of a highly satu- 
rated line with a small covering factor). Using 40% of 
the total flux for the main line, we flnd an HCO"'" col- 
umn density > 5 x 10^'*cm~^. This value is derived 
from the averaged monitoring data for 1996 and is 
slightly higher than the column density reported in 
Wiklind & Combes (1996a). Since the main line is 
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heavily saturated, the difference is duo to the much 
better signal-to-noise ratio of the new data. The col- 
umn densities have been derived using an excitation 
temperature Tx = 6K. Integrating the 3—3^2 spec- 
trum over the extent of the J=2<— 1 line sets a So- 
upper limit to the excitation temperature of 8K (a 
2<T limit implies Tx < 6K). 

3.3. The z=0.193 21 cm HI absorber 

We searched for redshifted C0(1-*— 0) absorption 
and/or emission as well as HC0+(2^1) absorption 
at z=0. 19267. Neither emission nor absorption was 
found. The CO and HCO+ spectra are shown in 
Fig. 4. 

3.3.1. Limits to the CO and HCO^ column densities 

The continuum level at the redshifted frequency of 

the CO(1<-0) and HC0+(2<-l) lines are 50 mK and 
26 mK, at the SEST telescope. However, when deriv- 
ing upper limits to the column density for molecular 
absorption we have to take the separate flux contri- 
butions into account. While the 21cm HI absorp- 
tion only occurs over the NE component (Lovell et 
al. 1996); putative molecular absorption can occur in 
front of both components. In Table 2 we give upper 
limits for three cases; (1) the total continuum flux, (2) 
60% of the flux, corresponding to the NE component, 
and (3) 40%, corresponding to the SW component. 
In all cases we assume an excitation temperature of 
6K. The ratio of the upper limits of the molecular 
column density to that of the HI is done for the NE 
component, since this line of sight causes the atomic 
absorption at z=0.193. 

We derived the HI column density from the spec- 
tra presented in Lovell et al. (1996), by fitting two 
gaussian profiles. The resulting column density is 
A/hi = 1-5 x lO^^TspZ/cm'^, where T^p is the spin 
temperature and / is the area filling factor. Since the 
HI absorption is only seen towards the NE component 
/ = 0.6/', where the factor 0.6 corresponds to the flux 
contribution of the NE component and /' is the area 
covering factor of atomic gas for this component. As- 
suming a spin temperature of 100 K, which is typical 
for atomic hydrogen gas in the Milky Way, the col- 
umn density of atomic gas is 2. -5 x lO^"/'"'^ cm^^. For 
the abundance ratios given in Table 1 we set f — 1. 

Using a typical value for 7V(H2)//^(CO) of 10*, the 
abundance ratio of molecular and atomic hydrogen is 
< 0.25 in the z=0.193 absorber. This limit is not very 



illuminating, since the column density of atomic gas 
is relatively low, corresponding to an extinction Ay « 
0.16 mag (e.g. Savage etal. 1977). Since CO becomes 
self-shielded only at Ay ^ 1 mag (e.g. Lucas & Liszt 
1994) we do not expect a signiflcant amount CO at 
these low HI column densities. Molecular hydrogen, 
however, can still be present, since it becomes self- 
shielded at lower A^hi- 

3.3.2. Limits to from CO emission 

Our CO(l-O) spectra of the z=0.193 HI absorption 
can also be used to derive an upper limit to the global 
molecular mass content of any intervening galaxy at 
this distance. The single dish telescope beam has a 
FWHP size of 52" at a frequency of 96.65 GHz, cor- 
responding to ~143kpc at z=0.193. Since atomic gas 
disks arc considerably smaller than this, an interven- 
ing galaxy causing the 21 cm HI absorption at z=0.193 
must be within our telescope beam. The effective 
beam of the interferometer data is 6" x2.4" (16.5x6.6 
kpc) with a noise rms of lOmJy/beam. The field of 
view is the same as for the single dish observation, 
since 15-m telescopes were used in both cases. 

The single dish CO spectrum shown in Fig. 4 has a 
channel-to-channel noise rms of 1.0 mK when binned 
to a resolution of 17kms^^. No emission is evident 
in the spectrum and we derive a 3a upper limit to 
the integrated CO intensity by assuming a linewidth 
of 200kms-i to be <0.17Kkms-^ The CO(l-O) 
spectrum does have a high frequency baseline curva- 
ture, which has been removed in the spectrum shown 
in Fig. 4. The uncertainty in the derived integrated 
intensity is limited by this curvature, but only for 
emission more extended than 1000 km s~^. Since such 
wide emission profiles are rarely encountered, we be- 
lieve that the baseline curvature is of little concern 
in this case. The interferometer data does not suf- 
fer from any baseline curvature and gives a 3(t upper 
limit to the integrated CO intensity of <0.09 K km 

The CO line luminosity, expressed in K kms~^ pc^, 
can be written as 

L'co = ^s*bDl '^^'^ 3 Kkms"^pc^ , (1) 

[1 + z) 

where Z>l is the luminosity distance, Iqo is the ve- 
locity integrated CO intensity expressed in Kkms~^, 
and il.s*b is the solid angle of the source brightness dis- 
tribution convolved with the telescope beam, Qf,. For 
distant sources and single dish telescopes, this convo- 
lution is dominated by the solid angle of the beam and 
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we can approximate it with Qs*b ~ i^b ~ 26.56 x B'^, 
(e.g. Solomon et al. 1997), where B is the full width 
at half power of the telescope beam measured in arc- 
secondfj^. The mass of molecular gas is obtained by 
assuming a conversion factor between CO luminosity 
and the column density of H2: = ct L'qq. We use 
a standard Milky Way conversion factor a = 4.6 Mq 
(Kkms-i pc^)"^ (e.g. Solomon et al. 1997). With 
these values our upper limits to /co transform into 
a 3cr upper limit to the H2 mass of 5 x 10^ and 
2.5 X 10^ M0 for the single dish and interferometer 
data, respectively. 

3.4. Additional absorption lines 

In addition to the 12 absorption lines presented in 
Wiklind & Combes (1996a) we have detected addi- 
tional absorption lines from CO(4^3), H2CO(2ii <— 
loi) and H^^CN(2^1). The spectra are presented in 
Fig. 5 together with a new spectra of H'^'^C0"'"(2<— 1) 
obtained with the IRAM 30-m telescope. 

For the first time we detect one of the CO lines 
in this system. The redshift of the absorber is such 
that the three first rotational transitions of CO fall 
at frequencies unfavorable for ground based observa- 
tions. The CO(4^3) line, however, is shifted into the 
1-mm band. We observed this line with the IRAM 
30-m telescope in May 1996 and the spectrum is pre- 
sented in Fig. 5a. The line is rather strong, with a 
depth almost equal to that of the saturated transi- 
tions of HCO+ and HCN. Using 40% of the total flux 
and an excitation temperature of 6K, the CO col- 
umn density is 2 x lO^^cm"^. Should be higher, 
say 20 K, the column density becomes 0.2 x lO^^cm"^. 
For higher excitation temperatures the column den- 
sity increases again. The main uncertainty in this 
derivation is the continuum flux. Due to unstable 
weather conditions it was not possible to derive T^. 
We have used the continuum flux determined from 
the 3- and 2-mm bands, which corresponds to 1.2 Jy 
and converted this to antenna temperature in the 1- 
mm band. Although indirect, the energy distribu- 
tion of the radio source PKS1830-211 is flat enough 
to make this approach viable. The excitation tem- 
perature of CO can be higher than the low values 
found for HCO+, HCN, N2H+, CS, etc. (see Wiklind 
& Combes 1996a). The reason is the lower electric 

■^Solomon et al. (1997) do not take the gaussian shape of the 
telescope beam into account and their molecular mass formula 
therefore gives a 13% lower mass estimate. 



dipole moment of CO, making it more easily collision- 
ally excited and giving CO a higher level population 
at rotational levels J>3 than do HCO+ and HCN. Dif- 
ferent excitation temperatures for CO and the other 
molecular species is consistent with the 'weak LTE' 
conditions assumed when deriving the excitation tem- 
perature, where it is assumed that the level popu- 
lation of each molecule is characterized by a unique 
temperature (that is: Tx = Tlot), but this tempera- 
ture can be different from one molecular species to 
another (see Wikhnd & Combes 1995, 1997). 

We also report the detection of ortho-H2CO through 
the JKaKb — 2ii ^ loi hne (Fig. 5b). The hne ap- 
pears to be non-saturated with a maximum optical 
depth ~1. The column density of ortho-H2CO is 
2 X 10^^ cm^^, assuming Tx = 6K and a flux frac- 
tion of 40%. 

Absorption of H^'^CO+ was reported already in 
Wiklind & Combes (1996a), but the present spectra 
represents a significant improvement to the signal-to- 
noise. The detection of this isotopic species was one of 
the main arguments in favor of a high H^^C0+(2^1) 
line opacity. The column density of H^'^CO"'", de- 
rived assuming the same excitation temperature as 
for the main isotope of 6 K, is 2 x 10^^ cm^^. The 
Hi2CO+/Hi3CO+ ratio is thus > 25, possibly much 
higher. From the similar absorption depth of HCO+ (2<— 1) 
and HCN(2^1) we inferred that also the HCN hne 
was heavily saturated (Wiklind & Combes 1996a). 
This is confirmed by the detection of the H^^CN(2^1) 
line with the IRAM 30-m telescope (Fig. 5d). The 
column density of H-'^^CN is 6x 10"'^^ cm~^ (again using 
40% of the total flux and adopting an excitation tem- 
perature of 6K). The Hi^CN/Hi^CN ratio is > 50. 

4. Discussion 

4.1. The lensing galaxy at z=0.886 

The detection of a second absorption line compo- 
nent at a redshift z=0.886, shifted in velocity relative 
to the main one and associated with the NE image of 
the background source, opens up a unique possibility 
to derive the dynamical mass of the central few kpc of 
the lensing galaxy. This is valuable information when 
modelling the lens, specifically when using time de- 
lay measurements to derive cosmographic parameters 
such as the Hubble constant. At the present this is 
hampered by the non-detection of the lensing galaxy 
at optical and infrared wavelengths (e.g. Djorgovski 
et al. 1992). This is likely to change with new and 
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more powerful infrared detectors. 

It is, however, possible to use very simple but pow- 
erful arguments, based on the assumption that the 
lensing galaxy is a 'normal' disk galaxy, to get rea- 
sonable ranges for the mass of the lensing galaxy as 
well as the redshift of the background radio source 
already from the present data. 

The observed line of sight velocity in a disk galaxy 
of inclination i, having a rotational velocity V{R) is 

Vohs = l^(-R) cos^sini , 

where R is the galactocentric radius in the plane of 
the galaxy and 9 is the angle from the line of nodes, 
as measured in the plane of the galaxy. R and 6 are 
related to the radius r and angle (p measured in the 
plane of the sky through 



R = r i cos 



2\ 1/2 



COSi^ 



, tan (p = tan 6 cos i . 



If we assume (it will be justified below) that i is small 
(i.e. i <J 30°), the differences between radii and an- 
gles measured in the plane of the sky and those in 
the plane of the galaxy are small (cosi ^ 0.87). This 
is a reasonable assumption since, as we will show be- 
low, the dynamical mass required to cause the im- 
age separation of PKS1830-211 implies that the rota- 
tional velocity of the lensing galaxy is more than twice 
larger than the observed velocity separation of the 
main and satellite absorption lines (e.g. sini < 0.5, 
and i < 30°). If we also make the assumption that 
the rotation curve is similar in shape to those found 
in nearby galaxies, we can assert that the rotational 
velocity is approximately constant at radii ^ 2kpc. 
This is seen from inspection of the rotation curves de- 
rived for all different types of spiral galaxies by Rubin 
et al. (1980, 1982, 1985). Two different lens models 
have been produced for PKS 1830-2 11 (Kochanek & 
Narayan 1992; Nair et al. 1993). They both agree 
that the center of the lensing galaxy is close to the 
center of the ring seen in radio interferometer data 
and that the lens has a small eccentricity. The model 
of Nair et al. places the center of the lens somewhat 
closer to the SW image than that of Kochanek & 
Narayan. In the former case the SW and NE cores are 
seen at projected distances of 1.8 kpc and 3.8 kpc from 
the center. In the model by Kochanek & Narayan, 
both cores are ~3 kpc from the center of the lens. In 
either case, it is reasonable to assume that the ro- 
tational velocity of the lensing galaxy has reached a 



relatively constant value at the projected location of 
both the SW and NE cores; V{R) « Vq. The observed 
velocity difference between two points in the galaxy 
is then 

AK)bs = Vo (cos^i — cos ^2) sini , 

where dj is the angle between the line of nodes and 
the line joining the center of the lens with the point 
where the line of sight crosses the disk, measured east 
of north. 

The angles 6, are 52° and 204° for the NE and SW 
components, respectively. This gives a simple relation 
between the maximum velocity Vq, the inclination i 
and the observed velocity separation of the absorption 
lines: \AV\ — Vo(1.5cosV' + 1.2 sin ■(/;) sini, where ip 
is the position angle of the major axis of the lens- 
ing galaxy. Geometrical considerations of the lensing 
morphology implies that the position angle is either 
~ 0° or ~ 90° (Kochanek & Narayan 1992; Nair et 
al. 1993). The mass inside a radius R is given by the 
virial theorem 



M(< R) 



AV^R 



7 X 10^ 

'■ 2 ■ 

sm z 



MfT 



I 



where (3 is the factor depending on the position angle 
of the lens: /3 = 1.5 for V = 0° and /3 = 1.2 for 
tp = 90°. We will adopt a position angle '0 = 0° in 
the following, keeping in mind that values can change 
by at most 20% should the position angle be 90°. 
In the last equality we used AV = 147kms~^ and 
R = 3 kpc. 

The two radio-components in PKS 1830-2 11 have 
almost identical substructure, suggesting that the 
components are two parity-reversed images of the 
background quasar caused by a nearly perfect align- 
ment of the lens and the source (e.g. Nair et al. 1993). 
The separation of the NE and SW images can there- 
fore to first order be described as the effect of a single 
point mass deflector of mass M, with the image sep- 
aration given by the 'Einstein angle' A^e 



A0T. 



GM 



5.7 X 10" 



arcseconds . 



In the last equality M is expressed in solar mass and 
the distance parameter D in Gpc; D = DdDg/Dds, 
where _Dd is the distance between the observer and 
the lens, Dg the distance between observer and the 
source, and D^s the distance between the lens and 
the source, all measured as angular size distances (e.g. 
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Schneider, Ehlcrs & Falco 1992). The last equation 
yields: 

Af w 3 X 10^" D Mq , 
[D in Gpc), and with the virial theorem: 

Vo ~ 220\/z? kms~i , (2) 

In Fig. 6 we plot as a function of the source red- 
shift for a fixed lens redshift of z=0. 88582. The dis- 
tances are derived using Ho — 75kms^^ Mpc^^ and 
for qo — 0.5 (full-drawn line) and qg — 0.05 (dotted 
line). D remains larger than 2 for all source redshift 
less than zsa7. This, together with the last equation, 
means simply that the observed image separation of 
~1" already requires a massive galaxy for the lens: 
for a redshift of the background source > 3, the re- 
quired rotational velocity is 350 km s^^, correspond- 
ing to a massive early type spiral. Lower redshifts 
for the background source would require higher rota- 
tional velocities, corresponding to an even more mas- 
sive galaxy, which is not realistic for spirals (unless 
several galaxies are considered to contribute to the 
lensing). Rotational velocities typical for late type 
spirals can be excluded, since the redshift of the back- 
ground source becomes forbiddingly large. Now, us- 
ing the relation AV' = 1.5Vb s'mi derived above, with 
AV = 147kms~^ and Vo = 350kms~^, yields an in- 
clination of i = 16°. 

How sensitively do these conclusions depend on 
the assumptions? (1) The approximation of the lens- 
ing mass as a singular point mass is only correct to 
first order if the lens and the background source are 
aligned. This is approximately true for PKS1830- 
211. (2) The derivation of inclination assumes that 
the absorbing gas is on circular orbits. If it has non- 
circular motion with an amplitude comparable to the 
rotational velocity, this simple approach is not valid. 
However, the lens model of Kochanek & Narayan 
(1993) independently finds a low inclination. (3) The 
presence of an absorber at z=0.193 introduces extra 
shear and, possibly, convergence of the ray bundles. 
However, given the relatively simple image configu- 
ration it is not likely that the extra mass concentra- 
tion at low redshift adds to the image multiplicity. It 
can, however, influence the interpretation of time de- 
lay measurements and Garrett et al. (1996) find from 
VLBI observations that the NE core component has a 
more extended and complex morphology than the SW 
component, which could be caused by the extra mass 
at z=0.193. Hence, our results regarding the mass of 



the lens at z=0.886 and the lower limit to the redshift 
of the source are quite robust and we can make the 
following postulates: 

1. the lensing galaxy is an early type spiral galaxy 
seen at a low inclination {i < 20°). 

2. The mass within a radius of '^3kpc is in the 
range ~ 6 — 9 x 10^° M©, which is similar, within 
a factor 2, to the Milky Way galaxy. 

3. The redshift of the background source is large, 
most likely Zg > 3. 

4.2. The covering factor of the SW compo- 
nent 

Frye et al. (1997) observed PKS1830-211 with the 
BIMA miUimeter array and found HC0+(2^1) ab- 
sorption only over the SW component, as we do. How- 
ever, the absorption line over this region did not reach 
zero level, indicating a filling factor of 70-80%. This 
would mean that the prospect of using single dish in- 
struments to monitor the continuum and depth of sat- 
urated absorption lines with the aim of deriving the 
time delay between the NE and SW components can 
lead to ambiguities. However, interferometric imag- 
ing of an object at low elevations is difficult in the 
millimeter band and it is unclear how the accuracy of 
the interferometer data can be assessed. The flux ra- 
tio between the NE and SW continuum components 
in the BIMA data is 1.14 ± 0.05| (Frye et al. 1997), 
compatible with our value at the same epoch, but 
their total flux is much higher that what we find, ei- 
ther with Plateau de Bure, IRAM-30m or SEST. In 
any case, the covering factor is quite high, almost 1 
for our measurements (see Table 1), and is also com- 
patible within 1 for the data of Frye et al (1997) when 
their noise level is taken into account or if their con- 
tinuum level is overestimated by a factor ~1.4. We 
do not think that this uncertainty is the main one in 
the time-delay determination, the major uncertainty 
being the continuum level measurement. 

It is not unexpected that a molecular cloud cov- 
ers entirely one of the PKS1830-211 image. The 
angular extent of the continuum emission region in 
the millimeter regime is not well known. The SW 
image remains unresolved in 15 GHz VLB A obser- 
vations (Garrett et al. 1996) at a resolution of 

^Taking the quoted fluxes and errors from Frye et al. (1997), 
we get a flux ratio of 1.14 ± 0.24. 
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0.6 X 0.2 milliarcscconds (mas). The NE component, 
however, is resolved with a complex structure. An 
angular extent of 0.2 mas corresponds to ~l.lpc but 
the size of the continuum source at millimeter wave- 
lengths can be considerably smaller. Millimeter VLBI 
of the BL Lac 3C446 found an upper limit to the an- 
gular extent of 30 /xarcseconds (Lcrner et al. 1993). 

In Fig. 7 we show the spectrum of the HC0+(2<— 1) 
absorption obtained with the Plateau de Bure in- 
terferometer with a velocity resolution of 0.5kms~^. 
The profile is asymmetric with an extension towards 
negative velocities. This part of the profile is either 
optically thin or has a covering factor less than the 
narrow component. Saturated absorption lines of a fi- 
nite width should be flat-bottomed when viewed with 
a high velocity resolution. The bottom of the narrow 
profile shown in Fig. 7 appears to be deeper at neg- 
ative velocities than at positive. The effect is only 
at the 1-2% level. However, a similar asymmetry 
is seen in the lower resolution HC0"'"(2<— 1) spectra 
obtained from both the IRAM 30-m and SEST tele- 
scopes (Fig. 2). This asymmetry could represent small 
differences in the covering factor of the saturated ab- 
sorption, corresponding to linear dimensions of -^5% 
of the extent of the continuum source. 

5. Summary 

We have reported further molecular lines in ab- 
sorption towards PKS1830-211, and in particular 
CO (4^3), the first CO line. Molecular abundances 
are compatible with Milky Way values. 

We have also carried out interferometric measure- 
ments at three different epochs, and found that the 
magnification ratio between the two images has de- 
creased from 1.8 to 1.0 during one year. These mea- 
surements confirm that the main V=0 absorbing com- 
ponent covers the SW image. We report about the 
detection of a second component at V=-147.5kms~^, 
which covers only the NE image. 

The presence of two absorption lines, one corre- 
sponding to the SW image of the background source 
and the other to the NE component, implies that the 
background radio source is situated at a redshift ^ 3, 
and is consistent with the lensing galaxy being an 
early type spiral seen almost face- on (i <; 1G°). A 
first determination of the time delay of 44±9 days has 
been reported by van Ommen et al. (1995). Given 
the model dependence of this measurement, it will be 
quite useful to try other independent determinations 



at various frequencies. We confirm here that with the 
SW component being nearly completely absorbed at 
3mm, single dish measurements are rapid and effec- 
tive for this (experiment. It will however be necessary 
to confirm from time to time that the filling factor 
is still 1 through interferometric measurement (but 
this has not varied in the past year). 

The lens geometry could be made more complex 
by the presence of another galaxy on the line of sight, 
detected by HI absorption by Lovell et al (1996) 
at z=0.1927. We detect neither molecular absorp- 
tion nor emission at this redshift; the galaxy might 
not be exactly aligned with the quasar, which would 
minimise its lensing effect. Alternatively, since the 
HI absorption is seen only towards the NE image, 
and molecular absorption has now been detected to- 
wards this component also, this z=0.1927 absorption 
could be re-interpreted as the absorption of a heavier 
molecule at z=0. 88489, the measured molecular red- 
shift towards this NE component. Exactly coinciding 
with this frequency we find the (7-1,0,7,6) transition 
of the PS molecule, in a low enough excitation state. 
This interpretation is not likely, however, due to the 
expected low abundance of PS in the ISM. 

We thank the staff at SEST and IRAM (both at 
the 30-m telescope and at the Plateau de Bure inter- 
ferometer) for efficient help and support during these 
observations. TW acknowledges support from NFR 
(the Swedish Natural Science Research Council) for 
this research and the European Southern Observatory 
for hospitality during a Visitor Grant Programme in 
Santiago. 
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Table 1 

Magnification ratios at different epochs 



Date Total Flux NE/SW Ratio /3, 



Oct 95 1.5 1.8 ± 0.1 1. 

Apr 96 1.45 1.2 ± 0.1 0.9 

Sep 96 0.93 1.0 ± 0.1 1. 



*The maximum covering factor of the SW com- 
ponent. 
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Table 2 

non detection of absorption at z =0. 19267 



Molecule 


Transition 


surest'' 




Tcont'^ 




Ntot [cm-2] 




N/NHi<i 






GHz 


GHz 


mK 


100% 


60% 


40% 




CO 


1+-0 


115.271204 


96.650 


50 


< 3 X 10^5 


< 6 X 10^5 


< 1 X 10i« 


< 2 X 10-5 


HCO+ 


2<-l 


178.375065 


149.559 


26 


< 3 X 10^2 


< 5 X 10^2 


< 8 X 10^2 


< 2 X 10-8 



*The rest-frequency of the observed molecules tciken from Lovas (1992). 

''Derived from i/obs = i^rost/(l + Za)- 
'^Continuum level measured in T^. 

''Ratio of molecular column density and HI (see text for details). 
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HC0+(2-l) 
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0.8 
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1 ' ' ' ' r 
HC0+(2-l) 



0.8 



0.6 



V (km/s) 

Fig. 1. — HC0+(2<— 1) spectra (bottom) obtained 
with the IRAM Plateau de Bure interferometer in 
September-October 1995 (dashed line) and April 
1996 (full line). The velocity resolution is 7.9kms~^ 
in both cases. Also shown is the relative location of 
the phase-center in RA and DEC (top and middle). 



0.8 



0.6 




Velocity (km s ) 



Fig. 2. — The HC0+(2^1) absorption profile from 
the 30-m IRAM (top) and the 15-m SEST (middle) 

telescopes, showing the second absorption component 
at — 147.5kms^^. The continuum level has been nor- 
malized to unity. The vertical dotted lines mark ve- 
locities at — 147.5kms~^ and Okms~^. In the bot- 
tom, we show the IICO"'"(3^2) spectrum, obtained 
with the SEST. No hue at — 147. 5km s^^ can be seen 
in the J=3- 2 line. The velocity resolution is 3.2, 2.2 
and 2.9 km s^^ from top to bottom. 
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o 

CD 



1.4 



1.2 



HC0(2-1) 
at z = 0. 88582 




09-apr 1996 



-500 





V (km/s) 



500 



Fig. 3. — Results of the fit of visibility data from April 
1996 observations: the integrated spectrum is at bot- 
tom, and the phase center position in right ascension 
at top: this shows that the V = — 147.5kms~^ com- 
ponent is absorbed in front of the NE image. 
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1.1 



0.9 
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0.9 



CO(l-O) PdB 



CO(l-O) SEST 




HC0+(2-l) SEST 




-300 



300 



Velocity (km s ^) 



Fig. 4.— The CO(l-O) and HC0+(2-l) spectra from 
the z=0. 19267 HI absorption system. No absorption, 
nor emission is present. The velocity resolution is 7.8 
and 10.5kms~^ for the CO spectra and ll.lkms"^ 
for HCO+. The velocity scale is heliocentric and the 
SEST spectra have been normalized to the contin- 
uum levels of 50 mK and 26 mK for CO and HCO+, 
respectively. 
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-200 200 -200 200 



Velocity (km s ^) Velocity (km s 

Fig. 5.- a)CO(4^3), b) H2CO(2n ^loi), c) 
Hi3C0+(2^1) and d)Hi3CN(2<-l) absorption at 
2=0.88582. All spectra except H2CO, have been ob- 
tained with the IRAM 30 m telescope. H2CO vi^as 
observed with the 15 m SEST. The velocity resolu- 
tion is 9.2knis~^ for the CO spectrum, 2.6knis~^ for 
the H2CO spectrum and 3.25 km s'^ for the Hi^CO+ 
and H^^CN spectra. The total continuum level has 
been normalized to unity and the velocity scale is he- 
liocentric in the restframe of the z=0. 88582 absorber. 
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2 4 6 

Redshift of source 



Fig. 6. — The distance parameter D = Dd-Dg/fds, 
where is the angular size distance between the 
observer and the deflector, the distance between 
the observer and the source, and £>ds is the distance 
between the deflector and the source. In this plot the 
redshift of the deflector is held fixed at z=0. 88582. 
The relation is given for a Hubble parameter Hq = 
75kms~^ Mpc~^ and for qo = 0.5 and qo = 0.05. 
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Fig. 7. High spectral resolution (0.5 km s"'^) 
Plateau de Bure absorption profile in HC0+(2^1) 
in April 96 (full line) and September-October 1995 
(dashed line). 
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